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INTRODUCTION 


In  reference  1,  page  21,  the  search  for  a  target  located  in 
the  plane  by  a  searcher  moving  in  space  is  studied.  The  instan¬ 
taneous  probability  of  detection  is  calculated.  It  isj 


y(R,H)« 


U) 


In  equation  1,  K  is  a  constant,  H  is  the  altitude  of  the 
searcher  and  R  is  the  distance  between  the  searcher' 3  projection 
and  the  target  (see  figure  1). 

The  constant  K  is  determined  experimentally.  According  to 
reference  1,  page  24,  the  sweep  width  for  the  detection  function 
(equation  1)  is: 


In  equation  2,  VR  is  the  speed  of  the  searcher  relative  to  the 
target.  The  sweep  width  W  can  be  measured  experimentally 
(reference  1).  Once  W  is  known,  K  is  found  by  inverting 
equation  2. 

This  inversion  is  not  easy  to  do.  Koopman  states,  "In  the 
majority  of  cases,  R  is  much  larger  than  H"  (reference  1,  page  9). 
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(3) 


i  - 

'  S  -i  __  ' 

^ 

He  then  replaces  equation  X  by: 

.  _  Y*  (R,H)  »  ^  ' 

'  —  -  •  ..  - 

and  equation  2  by  * 


The  integral  in  equation  4  is  done  by  a  change  of  variables.  Then 
(reference  1,  equation  29) 


In  order  to  obtain  equation  4,  it  is  assumed  that  H  «  X  .  But 
some  contribution  to  equation  2  comes  from  values  of  X  close  to  0. 
Also,  in  going  from  equation  2  to  4,  it  is  assumed  that  H  is 
small.  Aircraft  searches  are  typically  conducted  at  500  feet  to 
5,000  feet  altitudes.  Is  500  feet  close  to  zero?  Visual  searches 
from  ships  may  have  H  »  100  feet. 

In  this  paper,  the  integral  (equation  2)  is  evaluated  exactly- 
equation  5  is  then  approximately  correct  as  a  certain  dimensionless 
parameter  becomes  large.  Scaled  variables  are  used  to  introduce 
this  parameter.  Although  non-dimens ionalizat ion  is  familiar  in 
applied  mathematics,  it  is  not  used  much  in  operations  research. 

It  should  be, as  the  example  in  this  paper  demonstrates. 
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The  result  (equation  5)  is  correct  when  a  parameter  q  *  RH/VR 
is  large.  When  W/H  is  large,  q  is  large.  For  open  ocean 
searches,  W  is  typically  a  few  miles  and  H  is  typically  a 
fraction  of  a  mile.  Then  W/H  is  large.  There  are,  however, 
cases  in  which  W/H  is  not  large.  In  overland  search  or  search 
and  rescue,  when  the  search  is  conducted  over  a  forest,  W  may  be 
only  a  few  hundred  yards.  In  this  case,  W/H  may  be  small. 

In  the  second  section,  the  scaled  variables  are  introduced. 

The  third  section  relates  the  integral  (equation  2)  to  the  modified 
Bessel  function  of  the  first  kind.  The  validity  of  Koopman's 
result  is  then  verified  for  large  q.  In  the  fourth  section,  high 
altitude  or  small  sweep  width  searches  are  considered.  In  this 
case,  the  result  (equation  5)  is  not  valid.  Instead,  W  is 
proportional  to  H  .  Numerical  evaluation  of  W  is  given  in  the 
last  section. 


INTRODUCING  SCALED  VARIABLES 


Defina  ncp-dimens.ional  variables  w,  x,  q, 

W  *  Hw 

X  »  Hx  (6) 

q  -  K/HVr 

In  other  words,  height  H  is  used  as  the  length  scale  defining 
the  variables.  Equation  2  becomes 

+®  r 

w  *  J  1  -  exp  dx  .  (7) 

From  known  W,H  ,  we  can  find  w  .  From  the  exact  expression 
w  »  w (q)  derived  in  the  following  sections,  we  can  find  q  . 

From  q  ,  H  ,  VR  ,  we  have  K  . 
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*  CONVERSION  TO  A  BESSEL  Pit'.  M  f'\ ' . 
ASYMPTOTIC  BEHAVIOR 


The  transformation  x  =  tan  0  ,  dx  «  se<  't3  0  converts 


equation  7  to 

TT 

2 

/  2  \  * 

w(q)  =  / 

_JT 

1  -  exp  f-2q  cos^  0) 

—  ** 

A 


sec  8d  8 


(8) 


In  section  5,  tables  of  w(q)  as  a  function  o?  q  are  given. 
The  integrand  (equation  8)  is  even  with  respect  to  0  .  Thus, 

v/2 

r0 


jn/Z  2  \ 

w(q)  =  2J  (l  -  e"2q  cos  0j  soc20de 


(9) 


Differentiate  equation  9  with  respect  to  q  to  obtain 


dw  . 
dq  ~  % 


//2  2 
/  e~2q  COS  ed6  . 


(10) 


Notice  w(q)  is  a  strictly  increasing  function.  Using  the  double 
angle  formula 


cos2  0 


1  +  cos  2  8 


(11) 


Equation  10  becomes 


it/ 2 


P-  =  4e-1^ 
dq  „0 


f  e-  9  29d8  . 


(12) 
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Set  a  -  20  ,  so  that  <30  =  (l/2)da  .  Equation  12  becomes 

* 

dw  -  f  -q  cos  a, 

3q  =  2S  J0S  d“'  ■  (13) 

The  integral  in  equation  13  is  related  to  the  modified  Bessel 
function  of  order  zero,  Ig(q)  /  see  reference  2.  Thus, 


dw 

dq 


2Tre“q  Ig(q) 


(14) 


For  large  q  ,  IQ(q)  has  the  asymptotic  expansion  (reference  2, 
page  377) . 


i0(q) 


_ 9 

21  (dqi1 


+ 


225 

3 !  ( 8  q73 


Using  equation  14  in  equation  15  gives 


(15) 


dw  _  2n  (  1  9  225 

dq  /2^q  (  8q  2! (8q)5  3!83q3 


Integrating  equation  16  gives 


w  ~  /2tr 


(16) 


(17) 


Actually,  equation  17  contains  a  constant  of  integration;  we  show 
in  the  appendix  that  this  constant  is  zero. 
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Prom  equation  ;17 ,  we  h^V®  that  q-*-+«°  gives  large  w  ,  and  it 
is  shown  in  the  next  section  that  w*0  as  q-*-0  .  Consequently, 
when  W/H  =  w  becomes  large,  we  must  have  large  q  ,  so  equation  17 
holds.  In  the  original  variables,  for  W/H  large , 


W  ~  2 


- ? 

2ttV_H 

X\ 


The  result  in  reference  2  is  the  leading  term  in  the  expansion  of 

the  exact  solution.  It  is  not  necessarily  valid  when  the  altitude 

is  small.  Rather,  it  is  valid  when  ^  >>  1  .  It  is  not  possible 

W 

to  use  equation  5  to  determine  K  when  R  is  0(1)  or  smaller. 

To  determine  K  for  a  given  sweep  width,  find  the  value  of  q 
that  satisfies  equation  13.  From  this  value  of  q  ,  K  can  be 
found.  If  K  is  known  and  K/VRH  >>  1  ,  then  equation  5  gives 
the  sweep  width  W  .  In  the  fifth  section,  tables  of  w(q) 
are  given. 
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HIGH  ALTITUDE  OR  SMALL  SWEEP  WIDTH  SEARCHES 


For  a  high  altitude  search,  H  is  large  so  W/H  may  not  b$ 
large.  For  small  q  , 


dw 

<*£*»* 

dq 


2tro 


-q 


2  -q 

q  e  M 


,  ir  „4  ~q 
f  ^  q  e 


(19) 


e‘q^  .  (20) 

The  constant  of  integration  in  equation  20  is  zero  because 
equation  8  shows  w*0  as  q*0  .  Thus, 

/  -  K/HV  \ 

W  -  2tt11  (l  -  c  k  ) 

as  K/HVr  *>  0  . 


and  for  q-*-0 


w  ~ 


2  ir  ^1 


') 


NUMERICAL  EVALUATION  OF  W 


It  is  not  advisable  to  use  equation  5  to  determine  K  from  the 
sweep  Width.  Rather,  once  W  is  known,  K  should  be  determined 
from  equation  8.  Equation  8  can  be  integrated  numerically.  In 
table  1,  w(q)  is  given. 

A  more  exact  evaluation  of  w(q)  is  determined  by  power  series. 
Define  y  »  l/2ff  *  dw/dq  .  From  equation  18 


y  -  e~q  In(q)  . 


IQ(q)  satisfies  the  equation 


(21) 


2 


dX 


0  2, 


+  q  “d^  "  q  J0  "  °  ‘ 


(22) 


Then  y(q)  satisfies  the  equation 


q  +  (2q  +  1)  +  y 

dq  '  ! 


(23) 


Assume  that 


w(q) 
so  that 


E,  v“ 

n-1 


(24) 


Z  naci 

n-l  n 


n-1 


(25) 
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After  derivatives  are  evaluated  and  substituted  into 
equation  30,  coefficients  of  qn  are  set  equal  to  zero.  Then, 


At  q  **  10,  the  exact  and  asumptotic  results  differ  by  1  percent. 
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C'i  CO  <3-  LO  VO  f"  cC  C 


TABLE  1 


EXACT  AND  ASYMPTOTIC  VALUES  OF 
w(q)  «  W/H  • 


3_  EXACT  ASYMPTOTIC 


4,23 
6.58 
S  .29 
9,70 
10,92 
12.02 
13.02 
13,95 
14.82 
15.63 


5.01 
7.01 
8,68 
10.03 
11.21 
12.28 
13.26 
14,18 
15.04 
15.85 
22.42 
27.46 
31 .71 
35,45 


Overflow  problems. 


w(q) 

'  -i'  T 


7o  ERROR 

18 
.  7 


3 

2 

2 

2 

1 

1 
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APPENDIX  A 


EVALUATION  OF  THE  CONSTANT  OF  INTEGRATION  * 

! 

! 

To  see  that  the  constant  of  integration  in  equation  17  is  j 

zero ,  substitute  y  -  /q'cos  0  ,  dy  =  -/q  sin  0  d0  in  equation  9  | 

•'  -  i 

to  obtain  j 


(A-l) 


Then  notice  for  large  q  , 


0(1/ /q)  (A- 2) 


(A-3 ) 


and  no  constant  term  appears. 
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